High nitrogen-vacancy density diamonds for magnetometry applications 
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Nitrogen- vacancy (NV) centers in millimeter-scale diamond samples were produced by ir- 
radiation and subsequent annealing under varied conditions. The optical and spin relaxation 
properties of these samples were characterized using confocal microscopy, visible and infrared 
absorption, and optically detected magnetic resonance. The sample with the highest NV~ 
concentration, approximately 16 ppm (2.8 x 10^* cm~'^), was prepared with no observable 
traces of neutrally-charged vacancy defects. The effective transverse spin-relaxation time for 
this sample was T| = 118(48) ns, predominately limited by residual paramagnetic nitrogen 
which was determined to have a concentration of 49(7) ppm. Under ideal conditions, the 
shot-noise limited sensitivity is projected to be ~ 150 fT/\/Hz for a 100 /im-scale magnetome- 
ter based on this sample. Other samples with NV~ concentrations from .007 to 12 ppm and 
effective relaxation times ranging from 27 to over 291 ns were prepared and characterized. 

PACS numbers: (07.55.Ge) Magnetometers for magnetic field measurements; (61.72.jn) Color cen- 
ters; (76. 30. Mi) Color centers and other defects; (81.05.Uw) Carbon, diamond, graphite 

I. INTRODUCTION 

Optical magnetometers based on alkali vapor cells can measure magnetic fields with excep- 
tional sensitivity and without cryogenics; however spin-altering collisions limit the sensitivity of 
small sensors pQ. Paramagnetic impurities in diamond, on the other hand, are a promising sys- 
tem for millimeter- and micrometer-scale magnetometers, because they exhibit long spin coherence 
times (several ms j2]) over temperatures ranging from liquid helium to well above room temper- 
ature. Nitrogen-vacancy (NV) centers are particularly promising because they have a spin-triplet 
ground state and convenient optical transitions, allowing for efficient optical pumping and mag- 
netic sensitivity. Recently, single NV-centers were used for nanometer-scale magnetometry [SI 3] 
and magnetometry based on high-density ensembles of NV" centers was proposed [5]. Here we 
discuss preparation of diamond samples with parameters optimized for room-temperature, optical 
NV-ensemble magnetometers. 

The parameters for sample preparation most relevant to magnetometry are the concentration 
and homogeneity of the NV" centers, as well as the abundance of other impurities such as neutral 
nitrogen-vacancy centers (NV'') and substitutional nitrogen in the lattice. There is a wealth of 
techniques for growth, irradiation, and annealing of diamonds that can be considered (see, for 
example, Refs. [HI El El Ej and references therein). In this work, we prepared ten samples under 
different conditions and observed the optical properties at each stage of development using confocal 
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microscopy, visible and infrared absorption, and optically detected magnetic resonance (ODMR). 
The results form an initial matrix that will help identify the most important variables and pinpoint 
optimal ranges in the sample preparation process. 



II. EXPERIMENTAL PROCEDURE 
A. Sample parameters and irradiation 



# 


Company 


Synthesis 


[N] (ppm) 


Radiation 


Dose (cm ^) 


Anneal Temp. [lst,2nd,3rd] (°C) 


1 


Sumitomo 


HPHT 


< 100 


e~ 


4.0 X 10^8 


[700,875,1050] 


2 


Sumitomo 


HPHT 


< 100 


e~ 


9.8 X 10^^ 




3 


Element-6 


CVD 


< 1 


e~ 


4.0 X 10^^ 




4 


Element-6 


CVD 


< 1 


e~ 


8.0 X 10^^ 




5 


Element-6 


HPHT 


< 200 


e~ 


8.2 X 10^^ 




6 


Element-6 


HPHT 


< 200 


e~ 


9.8 X 10^^ 




7 


Sumitomo 


HPHT 


< 100 


e~ 


9.8 X 10^8 




8 


Element-6 


HPHT 


< 100 


e~ 


4.0 X 10^^ 




9 


Element-6 


HPHT 


< 100 


H+ 


1.0 X lO^'^ 


[500,875,1050] 


10 


Element-6 


HPHT 


< 100 


H+ 


1.0 X 10^^ 


[800,875,1050] 



TABLE L Sample characteristics. Initial nitrogen concentration, [N], is given as specified by manufacturer. 
Each annealing was carried out for 2 — 2.5 hours. 



Ten diamonds were fabricated commercially by Sumitomo and Element-6 using either chemical- 
vapor deposition (CVD) or high-pressure, high-temperature (HPHT) synthesis. The samples are 
millimeter-sized diamond crystals with an initial concentration of nitrogen impurities ranging from 
less than 1 ppm to 200 ppm (1 ppm corrsponds to 1.76 x 10^^ cm~^ in diamond). In the case 
of the HPHT samples, growth sectors with high and low nitrogen concentrations were identified 
and three of these sectors in Sample 5 were studied separately (see Sec. IIID|. The crystals 
were subjected to either proton or electron irradiation at varying doses. Electron irradiation was 
done at Argonne National Laboratory using 3.0-MeV electrons (e") from a Van de Graaff generator. 
Thermal contact with a water-cooled copper target was used to keep the temperature of the samples 
below 150°C during the irradiation process in order to reduce vacancy- carbon recombinations |10) . 
The irradiation doses were chosen so that the resulting vacancy concentrations would be within 

[nTB"3| 



The 



two orders of magnitude of the initial nitrogen concentration, as discussed in Sec. 
proton irradiation was performed at the Conditions Extremes et Materiaux: Haute Temperature et 
Irradiation (CEMHTI) laboratory using 2.4-MeV protons (H+), and the samples were maintained 
below 80°C during irradiation. Table [I] shows the specific conditions under which each sample was 
prepared. 



B. Optical characterization and annealing 

After the irradiation, optical characterization of each electron-irradiated sample was performed 
using confocal-microscopy at room temperature [H]. Light from a commercial solid-state 532-nm 
laser was spatially filtered using a single-mode fiber, coUimated, spectrally filtered, and combined 
into the collection path by a dichroic beamsplitter, exciting an optical transition associated with 
the NV" center's excited-state manifold [12] . The excitation beam was focused to a sub-micron 
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spot using a 0.6-numerical- aperture microscope objective with 40x magnification. The resulting 
photoluminescence (PL) was collected through the same objective, spectrally filtered to remove 
reflected laser light, imaged onto a pinhole, and analyzed. For spatial imaging, the microscope 
objective was mounted on a closed-loop triple-axis piezoelectric translation stage, and the collected 
photoluminescence was detected using a silicon avalanche photodiode, with filters used to select the 
NV" phonon sidebands (PSB) in the spectral range ~ 647-804 nm. For spectroscopy, the collected 
photoluminescence was sent to a spectrometer incorporating a cooled CCD detector, where the 
spectrum was recorded for various locations throughout each sample. The NV° and NV" zero- 
phonon lines (ZPL) were fit to Gaussian lineshapes after subtracting the local linear background. 
Imaging tests on single NV centers in an Electronic Grade ultra-pure CVD diamond purchased 
from Element-Six indicated an effective measurement volume of 5.2 /im^, and the single-center 
fluorescence spectrum provided the absolute scale for determining the NV" concentrations in the 
higher NV" density samples. 

Following irradiation and after the initial round of PL measurements, the eight electron- 
irradiated samples were subjected to thermal annealing at 700°C in the presence of a forming 
gas (0.96 atm Ar, 0.04 atm H2), in order to prevent surface oxidation, for 2 — 2.5 hours. The two 
proton-irradiated samples were annealed for two hours in vacuum, one sample at a temperature 
of 800°C (Sample 10) and the other at 500°C (Sample 9). Afterwards, each sample was analyzed 
with confocal spectroscopy at room temperature. This process was repeated for nine of the sam- 
ples at an annealing temperature of 875°C and then 1050°C in the presence of the same forming 
gas. After the final annealing, a Varian Gary 50 UV-VIS spectrometer and Nicolet Magna IR-750 
Fourier-transform infrared (FTIR) spectrometer were used to measure visible and infrared absorp- 
tion spectra for some of the electron-irradiated samples. Additionally, confocal microscopy was 
used for ODMR measurements and the effective transverse spin-relaxation times were determined. 

III. RESULTS 
A. Proton-irradiated samples 

Development of micro-magnetometers in certain geometries may benefit from the use of proton 
irradiation, where damage is confined to a smaller volume than with electron irradiation. In this 
section we present results of the optical characterization of two identical proton-irradiated samples 
(Samples 9 and 10). 

1. Stopping Range of Ions in Matter (SRIM) and spectroscopic measurements 

When charged particles traversing a medium come near to rest, the stopping power dramatically 
increases, causing a thin layer of radiation damage to develop. This effect is pronounced for non- 
relativistic ions in diamond. The ionization losses limit penetration to a thin layer near the surface, 
where radiation-induced vacancies in the carbon lattice are formed. In contrast, relativistic electrons 
make it through the 0.5 — 1 mm thick samples, and, consequently, the vacancies are distributed 
relatively evenly. 

Figure [T]^a) shows the results of proton-irradiation Monte-Garlo simulations made using the 
Stopping and Range of Ions in Matter (SRIM) software [13]. The diamond sample was modeled as 
a pure ^^G layer with density 3.52 g/cm^ and displacement energy 45 eV [HI [13 [TB]. For the H"*" 
energy of 2.4 MeV, the simulations predict a dense layer of vacancies, recoiled interstitial carbon 
atoms, and stopped H"*" ions at a depth of ~ 35 ^m, and almost no radiation damage deeper into 
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FIG. 1: (Color online) (a) Monte-carlo simulation of radiation damage as a function of penetration depth for 
2.4-MeV proton radiation, using the SRIM software, (b) Depth profile of NV~ centers, taken by normalizing 
the integrated ZPL intensity by that of a single NV center, after the first annealing (see text). The concen- 
trations for Sample 10 have been divided by a factor of 10 for visualization. Errors on the concentration 

[B2 (c) NV~:NV° ZPL intensity ratio as a function of depth after the 
The lightpower was 100 ^iW. For Sample 10, beyond a depth of ~ 36 /xm, the NV" ZPL 



estimates are discussed in Sec, 
first annealing 



became indistinguishable from the background. 



the sample. If vacancy-to-NV conversion is relatively even throughout the sample, then one might 
expect a similar depth profile of NV concentration. 

Figure [l](b) shows the estimated NV" concentration, taken by normalizing the integrated ZPL 
intensity by that of a single NV center, as a function of depth for the proton-irradiated samples after 
the first annealing. Figure [l|c) displays the NV":NV° ZPL intensity ratio under the same conditions 
as Fig. [T]^b) and is discussed later in Sec. |IIIA3[ In Fig. ^h), Sample 9 (after annealing at 500°C) 
shows a thin layer of NV" emission with a small peak at ~ 35 //m, followed by a rapid decrease 
in concentration with depth. The peak intensity, however, is not nearly as dramatic as in the 
vacancy profile predicted by SRIM. Sample 10 (after the first annealing at 800°C) shows a sharp, 
monotonic decrease throughout the first 35 /um that is not at all consistent with the vacancy profile 
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FIG. 2: (Color online) (a) Depth profiles of the phonon-sideband PL intensity for Sample 10 at a typical 
location (left) and at a sector boundary (right), (b) Deconvolved intensity as a function of depth for Samples 
1, 9, and 10 (SI, S2, etc.) after the first annealing at 700, 500, and 800 °C, respectively. The PL intensity 
is normalized by the incident laser power (5,153, and 1.9 nW for Samples 1,9, and 10, respectively). The 
signal for Sample 9 is multiplied by a factor of 100 for visualization, (c) PL deconvolution for Sample 9 after 
the first and third annealing. The PL intensity after the first annealing has been multiplied by a factor of 
30. The damped oscillations are artifacts of the signal processing caused by rapid changes in concentration. 



predicted by SRIM. Upon further measurements (see Fig. [2[a) in Sec. IIII A 2|), it was seen that 
the qualitative features in the depth profiles of Sample 10 varied greatly for different locations. 
It is worth mentioning that in this sample the NV~ concentration varied by more than an order 
of magnitude for different locations all at a depth of ~ 12 fim from the surface. There was one 
location that had a concentration nearly three orders of magnitude lower than the mean. As shown 
in the next section, these spatial inhomogeneities in NV~ concentration are related to the presence 
of multiple growth sectors, with widely varying nitrogen concentrations, in these HPHT-synthesized 
diamonds. 



2. Photoluminescence depth profiles 

The spectroscopic measurements were followed by spatially resolved PL measurements where 
the position of the focus was varied by translating the microscope objective with a closed-loop, 
triple-axis piezoelectric stage and the fluorescence, in the 647-804 nm range, was collected. Figure 
[2]ja) shows examples of raw data both for a typical scan (left side) and one around a sector boundary 
(right). After averaging over the lateral dimension, x, the resulting depth profile was deconvolved 
using a similarly-processed profile from a sample containing a very thin layer (^ 200 nm) of NV 
centers as the response function, in conjunction with an optimized Wiener filter [17j. Figure |2]^b) 
shows the results for three different samples: both of the proton-irradiated samples as well as 
Sample 1, a high nitrogen-content electron- irradiated sample, for comparison. These profiles were 
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taken after the final annealing at 1050°C and were post-selected to ensure that they were from a 
single growth-sector. As expected for relativistic-electron irradiation, Sample 1 has a relatively even 
depth profile throughout the sample. The gradual decrease in intensity with depth is consistent 



with absorption of the 532 nm laser light, as seen in Sec. IIIB4 Samples 9 and 10, on the other 



hand, exhibit emission only from a thin layer of NV~ centers near the surface. The edge of this 
layer is at approximately the same depth as the edge in the vacancy profile predicted by SRIM, 
~ 35 iim from the surface, but the sharp peak at the end of the layer is absent. Figure [2]^c) shows 
the depth profile in a single-sector region for Sample 9 after both the first and final annealing; the 
results are discussed below. 



3. Discussion of fluorescence depth profiles 



There could be several potential explanations for the deviation of the observed NV fluorescence 
depth profiles from the vacancy distribution predicted by SRIM: 

1. The ZPL intensity and the integrated PSB fluorescence are poor indicators of the actual NV" 
concentration. 

2. The vacancies diffuse tens of /im during annealing. 

3. The concentration of nitrogen is insufficient to convert most of the vacancies to NV~ centers 
at the stopping peak. 

4. Other products of the irradiation, either the interstitial carbon, stopped H"*" ions, or the 
vacancies themselves, inhibit NV" formation by forming stable defects with the vacancies. 



We now briefly discuss each of these possibifities and show that the first three explanations are 
not supported by our observations. 

1. In the He"*" ion implantation work of Ref. [2], a similar deviation from the SRIM vacancy 
profiles was seen in the NV" ZPL intensities for various radiation doses. The primary cause posited 
in that work was radiation absorption by vacancies and graphitization defects. In order to set a 
limit on the absorption by such defects in our samples, we monitored the intensity of the diamond's 
first-order Raman peak at 573 nm during the spectroscopic measurements whenever it was distin- 
guishable from the background. In every instance, the 573 nm peak intensity did not change by 
more than a factor of two over the first ~ 100 fim of the sample, and at least some of this attenu- 
ation was due to absorption by the NV" centers themselves (see Sec. IIIB4). To our knowledge, 
there is no defect present in large concentrations in these samples that absorbs strongly in the NV" 
fiuorescence region but not at 573 nm. The neutral vacancy center, GRl, for example, absorbs 
roughly equally in these two spectral regions [IHl HH]- Therefore a factor of approximately two 
is the maximum possible attenuation of the NV" fiuorescence due to absorption by other defects. 
The 35-/im peak height predicted by the SRIM vacancy depth profile is several orders of magni- 
tude larger than the peak heights of the NV" depth profiles in Figs. [l](b) and[2|b), so radiation 
absorption cannot explain this deviation. 

2. The results shown in Fig. ^c) give evidence that vacancy diffusion cannot explain the 
deviation from the SRIM vacancy depth profile. The magnitude of NV" fiuorescence in the depth 
profile for Sample 9 exhibits an increase of approximately a factor of 35 between measurements 
taken after the first and third annealing, signifying that vacancies are diffusing to form NV" centers 
during the second and/or third annealing. However, based on the shape of these depth profiles, 
and taking into account the artifacts introduced by the deconvolution, a limit of 6r ^ 1.5 ^m can 
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be placed on the spatial range of diffusion of vacancies in this sample during annealing at 875° C 
and 1050°C for a time of t = 2 hours each. This is not enough diffusion to account for the large 
difference in peak height predicted by the SRIM vacancy depth profile. Assuming the diffusion 
coefficient at 875°C is smaller than that at 1050° and that the vacancies are not annihilated by 
other processes during annealing, this corresponds to an upper bound on the diffusion coefficient 
at 1050°C of Dv ~ V(4t) < 4 x IQ-^^ cmVs. 

3. To address the possibility that insufficient nitrogen levels limit the overall NV" concentration 
at the stopping peak, a rough measurement of the levels of substitutional nitrogen available for NV 
formation was made by monitoring the ZPL of the neutral NV*^ center at 575 nm. Substitutional 
nitrogen is a donor in diamond, believed to be responsible for the negative charge on the NV" 
center, and consequently a shortage of nitrogen leads to a greater proportion of neutral NV'^ centers 
P [501 E]. Figure [l|c) shows the NV":NV° ZPL intensity ratio for Samples 9 and 10 for the 
same depths as in Fig. ^h). The high ratio for Sample 10 is inconsistent with donor depletion. 
The ratio for Sample 9 is lower compared to that of similar samples irradiated at lower doses 
(NV" : NV° « 1 — 40 (see Fig. |4](b)), but the absence of a sharp dip near 35 fim suggests that 
electron depletion cannot account for the deviation from the SRIM vacancy depth profile. 

4. The possibility that other radiation damage inhibits NV formation cannot be verified in 
this study, but it is not inconsistent with our observations. As seen in Figure [l](a) , the recoiled 
carbon atoms are the most abundant defect present in the diamond, with more than 30 times the 
nitrogen concentration. The deposited H"*" also comprise over five times the nitrogen concentration. 
At these levels, it is more likely that during annealing the diffusing vacancies will recombine, or 
form other color centers, with these interstitial atoms before they encounter a nitrogen atom. ESR 
measurements of the interstitial carbon color center, R2, suggest that some recombination occurs 
during irradiation unless the diamond is cooled well below room temperature [TUl [52] (not the 
case here), and that during annealing the remaining R2 centers become mobile at a threshold 
temperature in the 400 — 700°C range [TUl US 122] • Samples 9 and 10 are actually two halves of 
the same diamond, but their first annealing was done at different temperatures: 500°C and 800°C, 
respectively. The temperature of these samples was maintained below 80°C during irradiation, yet 
their depth profiles after these different annealing temperatures show similar deviation from SRIM 
vacancy profiles. This suggests that, if interstitial recombination is the mechanism responsible for 
this deviation, the onset of self-interstitial mobility occurs at or below an annealing temperature of 
500° C. This limit still falls within the 400 — 700° C range of values for the onset of self-interstitial 
mobility found in the literature. 

In addition to recombination with interstitial carbon atoms, vacancies are also known to combine 
with hydrogen atoms to form VH" [23] and NVH" Plj defects. As we were not setup to observe 
ESR where these paramagnetic centers are best detected, they were not studied here, but they may 
limit NV formation and dephase NV" spins in magnetometry applications. 

In summary, the proton-irradiation and annealing produced NV centers in a confined volume of 
depth ~ 35 //m. The depth profiles of the NV" center concentrations showed a significant deviation 
from the expected damage profiles. We considered possible causes of such a deviation and surmised 
that excess radiation damage likely inhibited NV formation. The resulting [NV"] distributions 
may actually be better for some micro-magnetometer configurations, but more work is required to 
understand the details of the radiation damage. 

B. Electron-irradiated samples 

In contrast to proton irradiation, irradiation with electrons leaves a relatively even distribu- 
tion of damage throughout the sample, enabling the preparation of millimeter-sized samples for 
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FIG. 3: (Color online) Fluorescence spectrum at each stage of annealing for Sample 2 (see Tab. |T]|, as 
well as the spectrum of a single NV center in the ultra-pure CVD diamond used for calibrations. The 
measurement volume is approximately 5.2 /im^ located ~ 12 /im from the surface. The important spectral 
features are labeled (as described in the text), (b) Lateral (xy) scans of the ultra-pure CVD diamond showing 
fluorescence in the 647 — 804 nm spectral region from two of the single NV centers used to determine the 
fluorescence collection volume. The range for all of the scans is 30 /zm x 30 iim 



ultra-sensitive magnetometry. Here we present and discuss the optical and magnetic-resonance 
characterization of electron-irradiated samples (Samples 1-8) prepared under varied conditions. 



1. Fluorescence spectra and calibrations 

Figure[3](a) shows the room-temperature 532-nm laser-induced fluorescence spectrum for Sample 
2, a 2 X 2 X 1 mm^ Sumitomo HPHT sample irradiated with electrons at a dose of 9.8 x 10^* cm~^ 
(see Table for each stage of annealing, as well as the spectrum for a single NV in a non- irradiated 
Electronic Grade ultra-pure CVD diamond. The collection volume was located near the surface of 
the sample in order to minimize the effects of laser and fluorescence absorption. Before annealing, 
the main spectral feature is a sharp peak at 573 nm which is the first-order Raman peak of the pure 
diamond lattice. Additionally, the broad second-order Raman feature is visible for the single-NV 
center in the 600 — 625 nm region. After annealing, the much brighter NV'^ and NV~ zero-phonon 
lines, at 575 nm and 638 nm, respectively, with linewidths of approximately 1.5 nm full width at half 
maximum, become prominent features in the spectrum. The monotonically increasing background 
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is due to the overlapping broad phonon sidebands. For more details on the spectral features, see, 
for example, Ref. [2S]- 

In order to quantify the NV concentrations, the fluorescence collection volume was determined 
by taking several lateral scans of the ultra-pure CVD diamond at different depths and recording 
the 647 — 804 nm fluorescence of single NVs, as illustrated in Fig. [Sj^b). A collection volume of 
~ 5.2 fim^ was calculated from these scans. This value, together with the integrated ZPL intensity 
at 638 nm from the spectrum of one of the single NV centers (Fig. [Sjja)), were used to determine 
the NV concentrations throughout this work. We note that concentrations obtained in this way 
are expected to be a lower bound, as the laser power used for these single NV scans was just 
above saturation, and the effects of PL quenching due to non-radiative energy transfer to nearby 
impurities p6j were not quantified. Nonetheless, as will be shown in Sec. IIIB4 and IIIB 5 this 
method produced concentrations values that are in close agreement with absorption-based methods 
[771 EE] • Furthermore, a fluorescence-based approach enables the measurement of NV concentrations 
in optically thin samples. 



2. NV concentrations 

Figurejij^a) shows the NV" concentration for all of the electron- irradiated samples after each stage 
of annealing. All concentrations were determined by normalizing the integrated ZPL intensity by 
that of a single NV center, and the collection volume was ~ 12 /jm from the surface. For the samples 
with high nitrogen content, the NV" concentration was already near, or at, its maximum value after 
the first annealing at 700 °C. The samples with < 1 ppm nitrogen did not reach a maximum until 
the 875 °C annealing. We note that since we do not completely distinguish between the effects of 
annealing time and temperature, it is possible that annealing the low nitrogen-content samples at 
700°C for twice as long would produce the same effect. 

In order to test reproducibility, the data were taken on two different days, using different lateral 
locations in the sample, and in between days the mirrors and sample mount had been adjusted. 
The samples were analyzed after the 700°C (but before the 875°C) annealing on both days. The 
discontinuity seen in Fig. [ij^a) between days shows that the uncertainty in determination of the NV" 
concentrations is about a factor of two, caused by the aforementioned conditions. We estimate that 
the total systematic uncertainty due to radiation absorption, uncertainty in the collection volume, 
etc. is less than the spread in concentration due to spatial inhomogeneity. 

The NV":NV° ZPL intensity ratio is plotted function of annealing temperature in Fig. 
gb). The ZPL intensity ratio can be converted to a concentration ratio by using the difference 



in the Huang- Rhys factors (Sec. IIIC) for the two charge states (NV":NV'' ZPL intensity ratio is 
approximately a factor of two smaller than the corresponding concentration ratio). The ratio does 
not appear to have a strong temperature dependence, but it varies strongly as a function of position 
in the sample and/or laser power, as evidenced by the large discontinuity between the two different 
days (corresponding to two different locations and laser powers). The incident laser powers are 



shown in the legend of Fig. |4](b) because, as will be discussed in Sec. Ill C changes in power lead 



to varying degrees of photo-ionization, which, depending on the individual sample characteristics, 
affects the ratio by an order of magnitude or more. 

The very high NV":NV'' ratio for Samples 1 and 2 provides evidence that many of the nitrogen 
atoms had not been converted to NV centers, since if this were the case, there would not be enough 
donors for making the negative charge state (Sec. Ill A 3). Conversely, the data in Fig. |4](b) show 



the effects of electron depletion on the low nitrogen-content samples, leading to a significantly lower 
NV":NV° ratio. Since the specific collection location and laser power between days were different, 
there is a large degree of uncertainty in the absolute values of the NV" :NV° ratio for each sample. 
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FIG. 4: (Color online) (a) NV~ concentrations of electron-irradiated samples as a function of annealing 
temperature, (b) NV~:NV'' ZPL intensity ratio of the same samples as a function of annealing temperature. 
The data were taken during two different runs on separate days. On each day the collection location for 
each sample was the same, but in between days the location changed. The laser power for the data on each 
day is indicated in the legend. Filled plot symbols are used for all samples with high nitrogen content and 
empty symbols represent the two samples with < 1 ppm nitrogen. Error bars are not shown here but are 
discussed in the text. 



However, for a given day and a given sample, the laser power and collection location were constant, 
so it can be concluded from the flat temperature dependence seen in Fig. [ijjb) for Day 2, and to a 
lesser extent also on Day 1, that the NViNV" ratio is not strongly affected by annealing. 

Figure [5] shows the NV" concentrations obtained after the final annealing at 1050°C as a function 
of radiation dose. As expected, the higher radiation dose generally leads to larger NV" concentra- 
tions. The initial nitrogen content is also an important factor, as evidenced by the factor of ~ 25 
difference between concentrations for Samples 3 and 8. These samples were subjected to the same 
radiation dose but their initial nitrogen concentration differed by approximately a factor of 100. 

In conclusion, the doses and annealing temperatures listed in Tab. [I] were sufficient to create 
samples with a wide range of NV" concentrations. Annealing at 700 °C was sufficient to obtain 
the maximum NV" yield for samples with high nitrogen content, while an additional annealing at 
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FIG. 5: NV concentrations near the surface as a function of dose for electron-irradiated samples after tiie 
final annealing at 1050°C. The data are the same as used in Fig. [4|a). 

875 °C was required for the low-nitrogen-content samples. The high-nitrogen samples exhibited 
a large NV":NV° ratio, while the low-nitrogen samples had a significantly lower ratio. Lastly, 
the concentration-measurement technique based on photoluminescence was reproducible to within 
approximately a factor of two even when different lateral positions in the sample were used. 



3. Nitrogen-to-NV conversion 



# 


[N] (ppm) 


post-irradiation [Vacancies] (ppm) 


final [NV"] (ppm) 


[NV-]:[N] 


1 


< 100 


45 


10 


> 10% 


2 


< 100 


110 


16 


> 16% 


3 


^ 1 


5 


.012 


> 1% 


4 


<1 


1 


.007 


> 1% 


5 


< 200 


90 


12 


>6% 


6 


< 200 


110 


8 


>4% 


8 


< 100 


5 


.3 


> 0.3% 



TABLE IL Sample Characteristics. Final NV~ concentrations for the electron-irradiated samples are tab- 
ulated along with initial nitrogen content (reported by the commercial providers) and projected vacancy 
concentrations. These values are used to set a lower bound on the nitrogen-to-NV~ conversion efficiency, as 
seen in the last column. The error in NV" concentration estimates, and consequently the error in conversion 
efficiency, is a approximately a factor of two, as discussed in the text. 

A figure of merit for sensitive high-density magnetometry is the conversion efficiency of sub- 
stitutional nitrogen to NV" centers, because any remaining paramagnetic nitrogen accelerates the 
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rate of spin decoherence (see Sec. |IIID4l). For irradiation temperatures around those used in 



this work (< 150°C), the rate of vacancy production for 2-MeV electrons throughout millimeter- 
thick samples was determined experimentally to be ~ 1.25 vacancies/e~ /cm in Ref. [TD] and 
~ 1.5 vacancies/e~ /cm in Ref. [22] and theoretically to be ~ 2.2 vacancies/e~ /cm |29| . Since 
the electrons are relativistic, the difference in vacancy production between 2- and 3-MeV electrons 
is small [2S], so we assume a vacancy production for the 3-MeV electrons used in this work of 
~ 2 vacancies/e"/cm. Using this value, along with irradiation doses from Tab. |l| Tab. [n| displays 
approximate initial concentrations of vacancies before annealing. The initial nitrogen concentra- 
tions reported by the manufacturers and the final NV~ concentrations from Fig. [s] are also shown 
and these values are used to set a lower bound on the nitrogen-to-NV~ conversion efficiency. Sam- 
ple 2, having converted ^ 16% of the initial nitrogen into NV~ centers, shows the most promising 
conditions for high conversion efficiency. For this sample, the initial nitrogen concentration was 
< 100 ppm and the irradiation dose was 9.8 x 10^® cm^^, corresponding to ~ 110 ppm vacancies. 
As discussed in Sec. IIIB4 it is likely this dose was insufficient to maximize NV" formation. 



4- Visible absorption 

Optical magnetometry with spin ensembles requires that the samples be reasonably transparent 
over the spectral range of interest. In diamond magnetometry based on NV" centers, zero-phonon 
optical transitions at 638 nm and 1046 nm [50], as well as their broad phonon sidebands, provide 
possible channels for detection. To measure the opacity in this range, a Gary 50 spectrophotometer 
was used to obtain transmission spectra for four of the electron-irradiated samples (Samples 1 — 4) 
after the final anneafing at 1050°C. The visible spectra for a subset of these samples were also 
measured after the first anneafing, and no significant differences in the spectra between these stages 
of annealing were observed. 

Figure [6] shows the room-temperature transmission curves for these samples. For all four samples, 
the transmission in the wavelength region ~ 700-1100 nm is approximately 66%. The loss of 
fight can be accounted for by Fresnel reflection; the index of refraction of diamond is 2.4, which 
corresponds to a reflection of 17% off of each diamond surface for light at normal incidence in air. 

In the wavelength region to the red of the ZPL and below ~ 700 nm. Samples 1 and 2 exhibit 
absorption that increases rapidly with decreasing wavelength. Absorption at these wavelengths is 
due to the significant population of NV" centers in the excited vibrational states at room tempera- 
ture. After subtraction of the overall background and absorption due to the ZPL (basefines in Fig. 
[6|, the ratio of the integrated absorption in this spectral region (638-700 nm) to the total integrated 
absorption is 11(2)% for Sample 1 and ~ 13(2)% for Sample 2. If these values are taken to be the 
Boltzmann-distributed occupancy of excited vibrational states at room temperature, they would 
give a rough estimate for the effective ground-state vibrational-energy splitting [31] of 56(4) and 
52(4) meV for the NV— center, respectively. These values are in tolerable agreement with previous 
measurements based on emission spectra, where the effective vibrational-spacing in the ground 
state of the NV" center was estimated to be ~ 60 meV [32] 155] . 

In the region around ~ 628 — 648 nm, inset in Fig. |6} the NV" ZPL appears prominently in 
the spectrum for Samples 1 and 2, where at the peak of absorption approximately 98% of the 
light is absorbed. This corresponds to four absorption lengths for this 1-mm path length, with 
approximately one absorption length attributed to just the ZPL and the other three due to overlap 
of the PSB. To the blue of the ZPL, at ~ 450 — 638 nm, resides a broad feature which is a 
combination of the NV*^ and NV" phonon sidebands. The ZPL for NV° at 575 nm is obscured by 
this strongly absorbing region. For Samples 3 and 4, there are no visible spectral features, due to 
the low NV concentrations in these samples (see Tab. |ll|. 
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FIG. 6: (Color online) Room temperature transmission curves for four electron-irradiated samples after the 
final annealing at 1050°C. The beam path (i.e. the sample thickness) was 1 mm for Samples 1 and 2 (SI 
and S2) and .5 mm for S3 and S4. The ZPL at ~ 638 nm is visible for SI and S2 and is blown up in the 
inset. The dotted lines represent the baselines used for Sample 1 in the effective vibrational-splitting and 
Huang- Rhys factor calculations (the baseline for Sample 2 was similar). 



Based on the absorption data, it is possible to estimate the electron-phonon coupling parameter, 
also known as the Huang- Rhys factor [31], S = —\n.{IzpL/ hot), where Izpl and hot refer to the 
integrated absorption intensities after normalization by the light frequency. While this relationship 
is only strictly correct at zero-temperature where quadratic electron-phonon coupling is suppressed, 
it was shown in Ref. [33] that, in NV" centers, the temperature-dependence of S is weak. In that 
work it was reported that S increases by less than 30% from 2 to 350 K. Analysis of the data 
in Fig. [6] gives, after subtracting the linear background (dashed lines in the figure), S « 3.9 for 
Sample 1 and 4.0 for Sample 2, which is within the range of values in the literature for NV" centers 
at or below 80 K (4.8,4.3,2.8 in Ref. [33], and 3.7 in Ref. [32]). The value used for hot is a 
slight over-estimate, as the phonon sideband from the NV*^ centers overlaps, but this effect is small 
(< 2%) since the concentration of NV° was small compared to that of NV" for both samples (see 
Sec. |IIIB2p . 

Another important observation from Fig. [6] is the absence of an absorption feature at 741 nm, 
the wavelength of the ZPL of the neutral vacancy center, GRl. The neutral vacancy concentration 
has previously been calibrated based on the GRl integrated ZPL in absorption [SHIES]- Based on 
this calibration and the lack of an identifiable ZPL, a conservative upper bound of 1 ppm can be 
placed on the concentration of neutral vacancy centers for Samples 1 and 2. This suggests that 
the NV conversion efficiency for these samples may be limited by the number of vacancies present 
after irradiation. Recall from Tab. [IT] that the vacancy concentration was predicted to be sHghtly 
larger than the nitrogen concentration for Sample 2. Since all of the neutral vacancies are gone 
after the final annealing, it is likely that there is a competing process to the vacancy capture by 
the substitutional nitrogen atoms. 

We also considered the possibility that many vacancies remain after annealing, but the vast 
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majority of remaining vacancies were actually negatively charged, due to the high concentration of 
nitrogen donors, and therefore do not contribute to the GRl absorption. The negatively charged 
vacancy center, NDl, is not visible in the absorption spectra because its ZPL is at 396 nm |25j. This 
wavelength coincides with the large background absorption, presumably due to photo-ionization of 
the substitutional nitrogen, prohibiting a direct measurement of the negatively charged vacancy 
concentration. However, as will be shown in Sec. Ill B 5 the concentration of positively charged 
nitrogen, the likely product of NDl formation, can be measured from infrared absorption spectra 
and was found to be 15(3) and 19(4) ppm after the final annealing for Samples 1 and 2, respectively. 
Since much of this concentration can be attributed to nitrogen's role as a donor in NV~ formation 
(recall that the NV~ concentrations were approximately 10 and 16 ppm for Samples 1 and 2, 
respectively), a conservative upper bound of 10 ppm can be placed on the negatively charged 
vacancy concentration in both samples. Since this upper bound represents a small fraction of the 
estimated initial vacancy concentration ([ll]), it remains likely that in order to maximize the NV 
conversion efficiency higher radiation doses should be explored. 

For Samples 1 and 2, the visible-absorption measurements can also be used as a cross-check on 
the NV" concentrations determined from the PL spectra (Sec. Ill B 2 ) [23 [2E] • Using the calibration 
pubHshed in Ref. [55], the integrated ZPLs for Samples 1 and 2 correspond to NV" concentrations 
of 8 and 10 ppm, respectively, which is in tolerable agreement with the concentrations measured 
from the PL spectra (10 and 16 ppm, respectively). As the calibration established in that work was 
intended for low-temperature absorption, the somewhat lower [NV"] values can be partly attributed 
to the slightly lower Huang- Rhys factor at room-temperature [5^ . 

In summary, visible-absorption measurements establish that even at high radiation doses the 
dominant source of absorption in the visible and near-infrared is due to NV centers, and not other 
radiation damage. The absence of GRl absorption suggests that higher radiation doses should be 
explored for optimizing NV yield. The absorption spectra were also used to confirm the NV" ground 
vibrational-state occupancy and Huang- Rhys factors from the literature and the NV" concentrations 
determined from the absorption ZPL provide good agreement with those determined from the PL 



spectra in Sec. IIIB 2 



5. Infrared absorption 

The presence of paramagnetic nitrogen has important implications for magnetometry, as dis- 



cussed below in Sec. |III D 4[ Infrared absorption measurements were made with a Nicolet Magna 
IR-750 FTIR spectrometer to measure the leftover nitrogen concentration after the final annealing 
(Fig. [?]). Nitrogen- rich diamonds display a broad vibrational absorption band peaking around 
1130 cm~^ [H^ that has been calibrated with respect to the substitutional nitrogen concentration 
using ESR and inert-gas fusion techniques [2ZIIM]- According to this cafibration, an absorption co- 
efficient of 1 cm~^ corresponds to a substitutional nitrogen concentration of approximately 22 ppm 
[55] . Figure [T] shows the infrared spectra for Samples 1 and 2 after the final annelaing, and a third 
Sumitomo HPHT diamond ("Yl") with very similar specifications that had not been irradiated or 
annealed. Based on the spectral feature at ~ 1130 cm~^, the nitrogen content for Samples 1, 2, and 
Yl were 58(8), 49(7), and 146(20) ppm, respectively, where the uncertainty is based on the spread 
of the cafibration values in the fiterature [ST] EH]- The sharp spectral feature at 1344 cm~^ present 
in all three samples is also associated with the single substitutional nitrogen [3S]. Both of these 
peaks are more intense for Sample Yl, presumably because Samples 1 and 2 contain NV centers 
and other nitrogen-related centers formed during irradiation and annealing. The most prominent 
feature in the FTIR spectra for Samples 1 and 2 is at 1450 cm^^, which was tentatively assigned 
to a vibrational mode of interstitial nitrogen [33 HO]- The smaller features at 1355 and 1358 cm~^ 
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FIG. 7: (Color online) Infrared spectra for Samples 1, 2, and "Yl" after the final annealing. Each spectrum 
represents averaging of 64 scans and the resolution of the spectrometer was set to 0.25 cm~^. The rapid 
oscillations are due to etaloning at normal incidence, and they have been smoothed by a moving average 
(black lines) in order to quantify the peak at ~ 1130 cm~^ (see text). Features mentioned in the text are 
labeled by their wavenumbers. 
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may also be related to this defect [H]. The production of this interstitial defect may be due to 
irradiation, since these features are not present in Sample Yl. 

The spectral feature at 1332 cm~^, present only in Samples 1 and 2, has been linked to positively 
charged substitutional nitrogen, N+ (Ref. [21]) • In that work, it was argued that because of 
nitrogen's role as a donor in NV~ formation, the concentration of N"*" after the final annealing 
should be roughly equal to the concentration of NV". Using the calibration of the 1332 cm~^ peak 
from that work (1 cm~^ of absorption corresponds to 5.5 it 1 ppm of N"*"), the NV~ concentrations 
for Samples 1 and 2 would be 15(3) and 19(4) ppm. Since nitrogen is a donor in the formation of 
other centers, such as the negatively charged vacancy (Sec. Ill B 4 1, it can be expected that [N"*"] 
would be somewhat higher than [NV"]. Thus these values for [NV"] are in excellent agreement with 
the values (10 and 16 ppm, respectively) measured using PL spectra in this work. 

In conclusion, the FTIR measurements indicate that a large concentration of paramagnetic 
nitrogen remains in the samples even after the highest irradiation doses and annealing temperatures 
studied here. The absorption of the 1332 cm~^ peak was found to be in excellent agreement with 
the NV" concentrations determined from PL spectra. 



C. Light-power dependence of NV*^ and NV" concentrations 

The neutral NV center, NV° presents a similar problem to ultra-sensitive magnetometry as 
the remaining substitutional nitrogen, since it is also paramagnetic [12] and therefore limits NV" 
dephasing times. Further, since the emission spectrum of NV*^ partly overlaps with the phonon 
sideband of the NV~ defect, the presence of NV° leads to a background fluorescence signal that 
reduces measurement contrast. In previous work [ISJ H^, it was shown that the NV" and NV° 
concentrations can change upon illumination by light in a range of wavelengths including 532 nm, 
the wavelength used for excitation in this work. In particular, in Ref. [33] it was reported that the 
NV":NV° ratio dropped with increasing laser power and identified direct photo-ionization of the 
NV" center and ionization of nearby potential nitrogen donors as the primary causes. 

Figure [s] shows the relative NV" and NV° concentrations determined from the PL spectra (Sec. 



IIIB 1) for (a) Sample 10 and (b) Sample 4 as a function of incident laser power. By relative, we 
mean that the NV" and NV'' concentrations have been normalized such that, at any given laser 
power, their sum is 1. In order to determine the NV'^ concentration, the difference in Huang- Rhys 
factors for NV*^ and NV" were used for normalization, since the NV'' ZPL for the single NV center 
(Fig. [3](a)) was not visible in the spectrum. As seen in Fig. [S] the laser-light intensity has a 
significant effect on the NV charge state for the high-nitrogen-concentration Sample 10, but does 
not affect the low-nitrogen-concentration Sample 4. Since at very low light power Sample 10 has a 
large [NV-]:[NV°] ratio and Sample 4 has a very low ratio, the results in Fig. [sjcan be interpreted 
as a small fraction of NV" being converted to NV'^ due to the laser light until the ratio reaches some 
minimum saturation level. Note that the range of laser powers explored for Sample 4 was higher 
than that of Sample 10, due to the low overall NV concentration in this sample. Similar plots 
were made after each annealing, but further annealing was found to have no measurable impact. 
While the mechanism for changing the NV charge state (i.e. one- or two- photon process) cannot 
be confirmed from this data set, the results indicate that sensitive magnetometry requires careful 
optimization of the laser light power. 
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FIG. 8: (Color online) Relative NV~ and NV" concentrations for (a) Sample 10 and (b) Sample 4 after the 
first annealing (700°C and 800°C, respectively) as a function of incident laser power. The concentrations 
were obtained from fluorescence spectra (insets) taken at a depth of approximately 36 fj,m, for Sample 10, 
and 12 /im, for Sample 4. Determination of NV° concentration was made by weighting the relative ZPL 
intensity according to the difference in room-temperature Huang-Rhys factors; S ~ 3.3 for NV° [25] and 
center (Sec. 



S - 4.0 for the NV 
0.5 fjxa for all measurements 



IIIB4I. The laser beam was focused to a diameter of approximately 



D. Optically detected magnetic resonance (ODMR) 

1. NV~ electronic structure 

One of the crucial features of nitrogen-vacancy ensembles for practical magnetometry is that 
the NV" centers can be optically pumped and interrogated using visible light over a broad range of 
wavelengths (see Fig. [6]). Figure [9|a) shows the crystallographic cell for the trigonal NV" center in 
one of four possible orientations, corresponding to the N-V axis along each of the dangling bonds. 
Figure ^h) shows the level diagram with energy levels and allowed transitions labeled according 
to the most current model posited in Refs. [SUl US]. Based on Cs^ symmetry |[22 and results 



18 




(a) (b) 

FIG. 9: (Color online) (a) Crystallographic cell for the NV center in diamond. The center has C^y symmetry 
and the line that contains both N and V defines the symmetry axis, of which there are four possible 
orientations. C-Carbon, N-Nitrogen, V- Vacancy, (b) Level diagram for NV~ center showing spin-triplet 
ground and excited states, as well as the singlet system responsible for intersystem crossing (see text). The 
transition wavelengths are for the ZPL, with radiative transitions indicated by solid arrows and non-radiative 
transitions by wavy arrows. 

from early electron-spin-resonance (ESR) and hole-burning experiments |17], it is known that 
the center's ground state is a spin-triplet and the state which can be directly accessed by an 
electric dipole transition is ^E. The ground state is split by 2.88 GHz due to the strong spin-spin 
interaction arising from the electron charge distribution being highly localized near each of the three 
carbon atoms At low temperature the state is an orbital doublet split by transverse strain, 
but at room temperature the orbital components appear to average [H], and thus it is represented 
by a single orbital with spin-spin splitting of 1.42 GHz |1S]. The ground state has a single 
non-degenerate electronic orbital and thus is immune to first-order shifts due to strain. 

In thermal equilibrium, the ground state sublevels are nearly equally populated, but after inter- 
action with light resonant with the — > ^E transition (excited here on the phonon sideband at 
532 nm), the \ms = ±1) sublevels become depopulated and the \ms = 0) sublevel eventually takes 
on the excess. Theoretical considerations and spin- dependent lifetime measurements suggest that 
the process responsible for the ground-state polarization is the decay, due to spin-orbit interaction, 
from the jm^ = ±1) sublevels of the ^E to a singlet state(s) with lower energy [121 SHI EO] . Pulsed- 
pump transient measurements have shown that the \ms = ±1) can decay back to the \ms = 0) 
through the singlet system after approximately 300 ns [12j. As this decay is mostly non-radiative, 
optical pumping into the \ms = 0) sublevel results in an increase in fluorescence with a contrast 
as high as 30%, limited by the branching ratio of approximately 0.5 for the \ms = ±1) sublevels 
decay to the singlet system [1]. Until very recently this singlet system, thought to be and/or 
^E from Cst, group theory considerations [121 EO], had not been directly detected. Weak infrared 
emission at 1046 nm was finally observed in Ref. [SU] and the magnetic- field and strain dependence 
of the fluorescence spectra indicate that the emission is indeed involved in the optical pumping 
mechanism. 

2. Optically detected magnetic resonance technique 



The transverse spin relaxation time, T|, is critically important for spin-precession magnetometry, 
as it places bounds on both the minimum sensitivity and laser power of broadband sensors (Sec. 
HIE). In order to measure T|, the technique of ODMR was employed after the final annealing. In 
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FIG. 10: (Color online) (a) ODMR spectrum for Sample 3. The 24 separate resonances correspond to two 
different |Atos| — 1 transitions, four different NV orientations, and three different ^''N hyperfine transitions. 
The orientations were split by a '--^ 50 G field, (b) Level diagram for the ground-state magnetic dipole 
transitions induced by the microwave field in the ODMR measurements. The selection rules for these 
transitions are lArn^l = 1 and Am/ = 0. The values for the hyperfine splitting are taken from Ref. 
|51) . (c) Microwave-power dependence (log-scale) of the relaxation rate for Sample 4 for a laser power of 
213 /iW. Microwaves from a signal generator were amplified and applied with a ~ 1-mm-diameter current 
loop with unmatched impedance. The microwave power was measured at the output of the signal generator 
prior to the ~ 30 dB amplifier. The FWHM linewidth of each |Atos| — 1, Am/ — transition for one of the 
orientations was plotted as a function of microwave power. A linear fit for the lowest-power points (dashed 
line) was used to extrapolate the linewidth to zero microwave power. The standard error in the fit parameter 
serves as the uncertainty in the measurement. 
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this technique, a microwave field is tuned to the ground-state ESR frequency of the NV" center and 
the fluorescence resulting from illumination by light exciting the transition is monitored 

[52] . Optical pumping transfers the NV" centers to the |ms = 0) ground state. An external 
magnetic field is applied, and the electronic spin sublevels \ms = ±1) couple to the component 
Bz along the N-V axis, leading to Zeeman energy shifts of AE = nisgsl-^BBz, where Qs = 2.0028 
is the NV" electronic Lande factor and ~ 1-40 MHz/G is the Bohr magneton. When the 
oscillating magnetic field has frequency equal to the splitting between the Im^ = 0) and \ms = ±1) 
magnetic sublevels, the component transverse to the N-V axis drives the magnetic dipole transition 
and repopulates the \ms = ±1) magnetic sublevels, counteracting the anisotropy produced by 
optical pumping. With the pump light still on, the increase in population in \ms = ±1) leads to 
non-radiative decay through the singlet system, reducing the fluorescence emitted by the center. 
Thus the ESR signal is the change in fluorescence registered when scanning the frequency of the 
microwave field. 

An example of an ODMR spectrum is shown in Fig. [lO|a). An additional layer of complexity 
in the spectrum is the presence of the hyperfine interaction with the nuclear-spin-one ^^N nucleus, 
which splits each lAm^l = 1 resonance into three hyperfine peaks (Am/ = 0) with an equal spacing 
of 2.3 MHz at zero-magnetic field [55]. Figure 10 b) shows the level diagram for the allowed ground- 
state transitions. The four different orientations of the NV" center are degenerate in the absence of 
magnetic field, but this degeneracy was removed with a static magnetic field of magnitude « 50 G, 
applied using a permanent magnet. Both fits to the data and numerical diagonalization of the 9x9 
spin hamiltonian ^I] confirmed that the spin-state mixing due to this applied field, particularly the 
component perpendicular to the N-V axes, causes the hyperfine splitting to vary by ^ 0.1 MHz. 



3. Results of ODMR measurements 



Sample # 


[N] (ppm) 


final [NV"] (ppm) 


Laser Power (/iW) 


T| (ns) 


1 


< 100 


10 


5 


66(8) 


2 


< 100 


16 


188 


118(48) 


3 


< 1 


.012 


130 


^ 291 


4 


<1 


.007 


213 


222(11) 


5 (spot 1) 


< 200 


1.0 


24 


27(4) 


5 (spot 2) 


< 200 


.02 


100 


128(64) 


5 (spot 3) 


< 200 


2.5 


240 


145(63) 


8 


< 100 


.3 


188 


114(25) 



TABLE IIL Effective transverse spin relaxation times, T2*, extrapolated to zero microwave power, for several 
electron-irradiated samples at the given laser powers (focused to ~ 0.5-/im diameter). For Sample 5, three 
locations, corresponding to spots in different growth-sectors, were investigated. The uncertainty estimates 
for all samples except Sample 3 come from the standard error of the linear fits, as shown in Fig. [Toj^c). For 
Sample 3, only one low microwave power, —30 dBm before the amplifier, was investigated (shown in Fig. 
[Toj^a)). The lower bound on T2* for this sample corresponds to one standard deviation below the statistical 
average of the resonance widths. 

The resonance profile corresponding to each N-V orientation and lAm^l = 1 coherence was fit 
with three Lorentzian functions of equal amplitude and equal linewidth, 7, given as full-width at 
half maximum (FWHM). Equal spacings between the centers of adjacent Lorentzians of 2.2 MHz 
were used, even when the ^^N hyperfine splitting could not be resolved. The effective transverse 
spin relaxation times, T| = l/(7r7), and associated uncertainties were obtained from the ODMR 
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spectra by extrapolating 7 to zero microwave power (Fig. [I0[c)) with a linear fit. Density-matrix 
calculations confirmed that 7 is, to lowest order, linear in the microwave power. The resulting 



r| times are tabulated in Tab. Ill for each sample at the given laser power. The laser-power 



dependence was not studied here though the light power was varied from sample to sample in order 
to maximize the signal on the photon counter. This dependence will be investigated in future work. 



4- Contributions to T2 

For the HPHT-synthesized Sample 5, ODMR spectra at three different locations, corresponding 
to sectors with very different initial nitrogen concentrations, were observed, and the NV~ concen- 



trations at these locations were estimated based on the integrated ZPL, as described in Sees. IIIB 1 



and IIIB 2 Since the irradiation was uniform throughout the sample, the wide variation of r|, 
from 27 ns to 145 ns is probably not due to vacancy-related defects other than NV centers, unless 
the formation of such centers is strongly correlated with the local concentration of nitrogen. 

For dilute spins, the contribution of NV"-NV" dipolar interactions to the magnetic-resonance 
broadening can be approximated by assuming that each NV~ center couples to only the nearest- 
neighboring NV" center. For an ensemble, this dipolar coupling leads to a spin-relaxation contribu- 
tion on the order of Tatv l/[(5s/^B)^?^Ary], where njw is the NV~ concentration [5]. For Sample 2, 
this corresponds to T^ry 1 jds, meaning that NV"-NV" coupling is significant but not necessarily 
the primary source of relaxation. For Sample 5, spot 3 has a much higher NV~ concentration but 
similar T| to spot 2, providing evidence that NV"-NV" relaxation is not the dominant source of 
relaxation in this sample. 

For nitrogen-to-NV" conversion efficiency of less than 50%, the contribution of substitutional- 
nitrogen dipolar interactions to NV" decoherence is expected to be even larger, as the paramagnetic 
nitrogen has similar coupling strength. Similarly to the estimate for NV~-NV~ interactions, the 
characteristic timescale for this decoherence is « l/{{gsHBYnN), where we ignore the small 



difference in g- factors between NV" and nitrogen |46j. Recall from Sec. IIIB 5 that Samples 1 and 
2 had final nitrogen concentrations on the order of 50 ppm, giving a value of T/v « 300 ns. Samples 
1, 2, and 5-8 all had similar initial concentration of nitrogen (< 100 ppm). Since the irradiation 
and annealing process converted no more than 16 ppm of nitrogen into NV~ centers, it is likely 
that relaxation due to remaining paramagnetic nitrogen is the dominant cause of relaxation in all 
of these samples. 

Another well-known cause for relaxation is inhomogenous broadening due to dipolar hyperfine 
coupling of nearby ^^C nuclear spins (natural 1.1% abundance in all samples). Based on experiments 
on single NV centers, this spin-relaxation contribution is in the range T^^^^ « 0.3 — 1.5 /xs [2][S5][M]. 
Those experimental values are not too far from the T| seen for Samples 3 and 4, where the nitrogen 
content was low enough that ^'^C is likely a major source of inhomogenous broadening due to dipolar 
coupling. 

In conclusion, dipolar broadening due to paramagnetic nitrogen is predicted to be the dominant 
dephasing mechanism in the high-nitrogen-concentration samples, followed closely by dephasing 
due to other NV centers and ^^C nuclear spins. For the low-nitrogen-concentration ([N]< 1 ppm) 
samples, dephasing due to the ^'^C nuclear spins was likely the dominant form of relaxation. 
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E. Prospects for ultra-sensitive magnet ometry 

The performance of a magnetometer based on electronic spins is, in principle, limited by the 
spin-projection noise. The minimum detectable magnetic field, 6B, is [Ij: 



where R is the measurement contrast, rj is the detection efHciency, N is the number of spin centers, 
and t is the integration time. It is interesting to estimate the magnetic sensitivity of a sensor based 
under ideal conditions. We choose Sample 2, as it has the highest product of density and T|, and 
we speculate that improved experimental techniques will allow the signal from nearly all NV centers 
to be collected with contrast R ~ 0.3 and efficiency t] ~ 0.1. Using these values the quantum shot- 
noise limited sensitivity would be ~ 150 fT/VHz for a 100 /im-scale magnetometer based on this 
sample. This approaches the sensitivity of superconducting quantum interference devices (SQUIDs) 
with similar spatial resolution [SS]. Moreover, the isotropy of the system makes it easily size 
scalable, and the advantage of a such a compact device working in a broad range of temperatures- 
from room temperature down to liquid helium temperatures-makes NV-ensembles an attractive 
candidate for many types of magnetic-sensing applications. Of course, much work needs to be done 
to experimentally realize the shot-noise limit, including thorough analysis and mitigation of sources 
of decoherence and improved experimental techniques. In particular, we note that the optimal 
optical pumping rate for this model magnetometer is on the order olNjqy /T2 = 2.8 x 10^^/(118 ns). 
This corresponds to approximately 8 W for 532 nm laser fight. Thus, lengthening the effective 
relaxation time is essential in order to improve sensitivity while maintaining practical light powers. 



IV. CONCLUSION 



Proposals for high-density ensemble magnetometers predict sensitivities rivaling the most sensi- 
tive atomic and superconducting magnetometers [5]. Development of such magnetometers requires 
reliable production of diamond samples with control over the NV" concentration as well as other 
defects which may contribute to spin relaxation. In this work, ten diamond samples were prepared 
under very different conditions and their optical and spin-relaxation properties were character- 
ized. The compiled results show that electron- irradiated samples annealed at 700° C for two hours 
can have NV~ concentrations on the order of 10 ppm, a reasonable benchmark for 100 /xm-scale 
magnetometry at the 100 fT/\/Hz sensitivity level. For samples with large initial nitrogen concen- 
trations approaching 100 ppm, the highest irradiation dose used in this work, ~ 10^^ e~/cm^, was 
still insufficient for maximum nitrogen-to-NV~ conversion. Nitrogen concentrations on the order 
of 50 ppm remain in the samples, so in order to maximize the conversion efficiency, higher doses 
will be explored in future work. Conveniently, the irradiated samples remained mostly transpar- 
ent in the spectral region of interest-absorption in the visible and near-infrared spectrum was due 
almost entirely to NV centers. Effective relaxation times above 100 ns have been observed for 
some of the high-NV~ concentration samples, likely limited by remaining paramagnetic nitrogen. A 
100 /xm-sized magnetic sensor based on this sample operating under ideal conditions should achieve 
a sensitivity of ~ 150 fT/\/Hz. In order to improve on this limit, future work is necessary to 
mitigate the dominant causes of relaxation. 
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